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We investigate whether baryogenesis is possible at a second order elec- 
troweak phase transition. We find that under rather general conditions, the 
departure from thermal equilibrium is suppressed by the expansion rate of 
the Universe, and hence baryon production is also suppressed by the ex- 
pansion rate. We conclude that if no additional sources of departure from 
thermal equilibrium such as topological defects are present, then electroweak 
baryogenesis is ruled out if the phase transition is second order or crossover. 
However, a non-vanishing net baryon to entropy ratio is generated, and we 
provide both upper and lower bounds for the result. Our technique is also 
applicable to other baryogenesis mechanisms taking place during and imme- 
diately after a second order phase transition. We estimate the lowest value 
of the transition energy scale for which the resulting baryon to entropy ratio 
might be large enough to explain observations. 
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1 Introduction 



Recently the constraints on the Higgs sector, at least in the Minimal Standard 
Model have improved so that the prevalent view today is that the electroweak 
phase transition is either very weakly first order, second order, or perhaps 
even crossover. The first evidence is a direct lower limit for the mass of 
the minimal Higgs particle > 65GeV (95% CL.) set in experiments at 
CERN (for a review, see for example 0) and based on the lack of Higgs 
production in hadronic decays of the Z boson. An indirect limit arises from 
the precision electroweak data; when the CDF and DO measurements of the 
top quark mass are taken account of (m^ ~ 181 ± 12GeV), one arrives at 
the la allowed range for the Minimal Standard Model (MSM) Higgs mass: 
26GeV< Mh < 230GeV (68% C. L.) §. In addition the requirement for 
(meta)stability of the electroweak vacuum p, implies a lower bound of 
Mh > 116GeV, and the requirement that the Standard Model couplings 
remain perturbative up to a scale Ap ~ Mpi ~ lO^^GeV results in a pertur- 
bative upper bound < 190GeV 0. Taking into account all of the above 
we see that there are strong indications that the MSM Higgs mass is between 
100 and 200 GeV. 

In theories with an extended Higgs sector (with two or more Higgs dou- 
blets) the experimental limits on the mass of the lightest Higgs scalar are 
somewhat weaker. For example, in the Minimal Supersymmetric Standard 
Model (MSSM) the direct limits from the LEP experiments on the light- 
est neutral Higgs boson mass are roughly Mh > 40GeV. The direct con- 
straints on the mass of the second neutral Higgs boson M^ become rather 
weak once the radiative corrections from other supersymmetric particles of 
unknown mass are taken into account. The meta-stability bounds do not 
apply to the MSSM vacuum, while the intrinsic upper bound on Mh is about 
140GeV (for < IQOGeV) §. So, in the case of the MSSM, the bounds on 
the lightest Higgs mass indicate a somewhat lighter Higgs particle: 40GeV 
<Mh< 140GeV. 

Next we discuss how the nature of the electroweak phase transition de- 
pend on the details of the Higgs sector, and in particular on the mass of 
the (lightest) Higgs particle. We now understand the nature of the phase 
transition in the Minimal Standard Model for a moderately light Higgs 
Mh < TTLw — 81 GeV. The pictures of the phase transition emerging from 
two- loop calculations [l^ , |ll[ and lattice calculations [|12| agree and indicate 
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that the transition is strongly first order for small Higgs mass and becomes 
weaker as the mass becomes comparable to mw- It is believed that for a 
sufficiently large Higgs mass the phase transition eventually becomes second 
order [|13[, [1^, or even a crossover, since at large Higgs masses the the- 
ory simply resembles more and more scalar theory. At what Higgs mass this 
exactly happens is not clear at the moment. Also how the dynamics of the 
phase transition is affected by the presence of more than one Higgs doublet, 
as is the case in the super symmetric versions of the Standard Model, is not 
completely understood. It is possible, for example, that the phase transition 



proceeds in two stages [0 , but we will not consider this possibility here. We 
are mainly interested in how the strength of the phase transition changes; 
even though there is no complete analysis of a general supersymmetric two 
Higgs doublet model, in the case of the MSSM it is known [jl^ that the stop 
sector tends to strengthen and a light CP odd neutral Higgs scalar tends to 
weaken the phase transition. We conclude that, given the above consider- 
ations, it is likely that the electroweak phase transition is second order or 
crossover. 

In this letter we will address the question: Is baryogenesis at a second 
order electroweak phase transition ruled out? (For reviews of electroweak 
baryogenesis see [|18|, [l^.) We will assume that no other sources of departure 
from thermal equilibrium are present except for those caused by the phase 
transition and the expansion rate of the Universe. This means, for example, 
that we will not consider the case when one has a network of cosmological 
defects |^ floating around that could drive the system out of equilibrium. 

The letter is organized as follows. In Section 2 we derive a simple dissi- 
pative equation of motion which out of equilibrium may lead to a net baryon 
number production. In the subsequent section we analyze this equation, de- 
rive an upper bound on the number of baryons produced and find that it is 
suppressed as the expansion rate of the Universe. In Section 4 we discuss a 
lower bound on the strength of baryogenesis at a second order phase tran- 
sition and use the result to estimate the minimal energy scale for a second 
order baryogenesis mechanism without extra out of equilibrium effects such 
as topological defects. In the final section we summarize our results and their 
consequences for electroweak baryogenesis. 
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2 Equation of motion 



We start from a standard near equilibrium formula of statistical mechanics 
for the rate of approach to equilibrium of some charge Q 

Q = -Tq^ (1) 

where ixq = AF, the chemical potential for Q, measures by how much the 
free energy changes when Q changes by one unit, and Tq is the rate of decay 
of Q (per unit time). This equation is a macroscopic form of the charge 
(non) conservation and can be thought of as an integral of the Boltzmann 
equation (see for example [Q). 



When applied to the baryon number {Q = B), we get (see |^T]] for a 
similar analysis): 

nB^^ = -^^^^^AB{sph), AF^^s=i)=fiB (2) 

B is the total rate of change of baryon number, Tsph/V the sphaleron rate 
per unit volume, T is the temperature of the plasma, /is = AF(ab=i) the 
chemical potential for baryon number (which we defined as the change in 
free energy when baryon number changes by one unit), and AB{sph) is the 
change in baryon number per sphaleron transition. In the symmetric phase 
the sphaleron rate is Tgph/V = /€(a^T)^ with k, ^ 1.1 |2^, |2^, and in the 
broken phase 

^ocexp-^, (3) 

where Esph is the sphaleron energy. 

We will now illustrate how using Eq. (|^) and a one loop formula for 
the free energy of fermions in a thermal plasma one can obtain an equation 
for near equilibrium baryon number production. The (one-loop) thermal 
contribution to the free energy reads: 

^ = ^E/ ^ln(l + exp(-/?(E,-/i,))), (4) 

where the sum J2i is over all fermionic degrees of freedom and the energy in 
the presence of a non-vanishing 6 is given (in the WKB limit) by 



mTgeef + m^Y^' for = ±1 (5) 
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where gg = — f|/2(f^ + f|) is the '^-charge,' (which is related to the axial 
charge), vi and V2 are the vacuum expectation values (vevs) of the two Higgs 
fields, respectively, and is the spin of the particle. In the zero mass limit 
the particles with = +1 reduce to the left handed fermions or left handed 
anti-fermions (anti-particles of the right handed fermions), whilst those with 

= — 1 reduce to the right handed fermions or right handed anti-fermions. 
Even though this relation is derived using equilibrium techniques and 6 is 
a time dependent quantity, we expect it to be a good approximation to 
description of near equilibrium processes since the typical time scale at which 
6 changes is given by the expansion rate of the Universe and the processes 
which equilibrate the system are typically much faster. 

For any particle species denoted by a subscript i, it follows from Eqs. ([^) 
and (0) that the number density rii obeys the following equation (c/. p6|): 



= -^^^iY^^jf^j (6) 

where uj are the stechiometric coefficients of the reaction and fij are the 
corresponding chemical potentials. In the case of sphalerons, Tsph/2V is 
the rate per unit volume for each of the following two sphaleron transitions: 
tLtibLTL--- ^ 0, tLbibi^T--- ^ 0, where dots (...) denote the particles of the 
lighter two families. The coefficients Ui of the first process are for example 
Utj^ = 2, i/ft^ = 1, etc, so that (|) for the left-handed top quarks and left- 
handed bottom quarks yields: 

hb^ = -^^[4:fXtL + 5/i6L + f^TL + (8) 
+4/ici, + + fJ'f^L + 2Ati/^ + 4/i„^ + 5fidL + f^cL + 2/ii.e] 
Having in mind that the total baryon number density B is 

B = l E (9) 

quark species 

and "hi^ = for right handed particles, one gets: 

B = -Np^^[3fitL + S/^bi + ^^TL + ^J'l^r (10) 
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where Np = 3 is the number of famihes. Note that we have ignored the 
anti-particles; it is rather trivial to include them in Eq. ([Til): one should just 
subtract the chemical potentials for the anti-particles. 

The final step is to relate the chemical potentials in this relation to the 
number densities. The particle number density can be easily obtained from 
Eq. d) as follows: 

n. = -^F=f^ - (11) 

dE, J (27r)3l + exp[/5(Ei-/i,)] ^ ^ 

Next we expand the number density to linear order in /Xj: 



After some algebra (see Appendix A) we get (to leading order in 9 and in 
the high temperature limit T > rrii) 

/i. = ^[{n, - n°) ± cim^)goe] , c(m?) = ^m,^(l - ^). (13) 

We have displayed the cubic mass term just to get a feeling when the high 
temperature expansion breaks down; for mj ~ T we expect it still to be 
reasonably accurate. Eq. (|lT]) can now be re-written as 



i 

where Np = 3 is the number of families, J2i ^iL = '^l is the total left handed 
fermion number density, which can be recast in terms of baryon and lep- 
ton numbers as "iBp + Ll. (In the above the equilibrium contribution from 
particles has been cancelled by the contribution from anti-particles 

— X^j^i*-) We write the final form of the baryon number equilibration equa- 
tion as follows: 

B = -QNp^ [3Bl + Ll + 6cim^,)gee] (15) 
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where the contribution from the hght fermions to c{mf) has been ignored in 
comparison to the top quark. The coefficient 6 = 2Nc in front of the 6 term 
is the number of degrees of freedom for the left-handed top quark. 

In the next section we discuss the solution to this equation assuming that 
the phase transition is second order or crossover. 

3 An upper bound on the baryon to entropy 
ratio 

We will now analyze Eq. (^) and, assuming that before the transition the 
Universe is in thermal equilibrium so that baryon number is zero initially, we 
will investigate what is the maximum baryon number produced at a second 
order phase transition in a two Higgs doublet model. 

First we simplify Eq. ([T5|) . Since the time scale for the weak sphaleron 
transitions is at least Tgph ~ 1/ {a'^^T) , r^ph is much larger than the time scale 
of strong sphaleron transitions Tgs ~ l/(a^T), so that at any moment when 
discussing only the dynamics of the weak sphaleron processes, the strong 
sphalerons will be to a very good approximation in chemical equilibrium, 
which means Bi = Br- We also have to make an assumption about the time 
scale for equilibration of the left and right lepton numbers tlr ~ 1/ {awnfT) 
(here yi ^ t/r ^ 10~^ is the Yukawa coupling constant). For simplicity 
we will assume t^r << r^p/j, which imphes that to a good approximation 
Ll = Lff. This will certainly be overwhelmingly satisfied when the weak 
sphaleron rate is exponentially suppressed, the case of most interest to us. 
When the symmetry is restored and when the sphaleron rate is un-suppressed, 
tlr > Tsph may be a more reasonable approximation. We will not study this 
case here since, as we will see below, most of the baryons are produced at 
a time when the sphaleron rate is exponentially suppressed. In addition we 
know that B — L is conserved by the the Standard Model. For definiteness we 
will set B — L to zero, which is motivated by the symmetric initial conditions. 
To summarize we have: 

Bl = Br Tsph » Tss (16) 

B = L (symmetric initial condition) 

Ll = Lji Tgph >> Tlr 
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/^From these constraints we easily infer Bl = Ll = B/2 which imphes 3Bl + 
Ll = 21?. Eq. (|15D can be now written in a simple form 

B = -TiB + Ce), C = ge^^{l~^), T = 12iV^^, (17) 

where rrit = ytip is the top mass. Here, is the Higgs expectation value which 
changes during the electroweak phase transition from = to = 0(T = 0). 

The solution to this equation with the zero initial baryon number can be 
easily obtained by the Greens function method 

B{t) = - f dt' Ce{t')T{t') exp[- fv{t")dt"]. (18) 

tin i 

In order to derive an upper bound on -B(t), we write 

Bit) < Bo f dt'T{t') exp [- /* dt"T{t")] = Bq (19) 

J tin 

with 

So = max,,,[,^„,,]|C(t')^(Ol (20) 

and we have assumed that /j*^ Tdt >> 1. This formula is our upper bound for 
baryon number production. Next we will show that, in the context of a simple 
model for the effective potential near the phase transition, Bq is suppressed 
by the expansion rate of the Universe. We believe that our argument is not 
limited to the simple form for the potential we use: 

V{^,T)=D{T'-T^)<P' + ^<P' (21) 

Here we quote the one-loop Standard Model values for the parameters: Tq ~ 
1.39M^ + (60GeV^)^ is the spinodal temperature, D ~ 0.18, and A depends 
on the Higgs mass (M^/2wq, vq = 246GeV) and, in addition, it is weakly 
(logarithmically) dependent on temperature. A ~ 0.05 — 0.16, depending on 
the Higgs mass, the lower bound corresponding to Mh = 80GeV, the upper 
to Mfj = 140GeV [^. The value of A is weakly dependent on temperature; 
for example when temperature decreases by about 10% (which will turn out 
to be the value at which (^) is maximized), the value of A changes by not 
more than 10%, which we can ignore. Since we do not trust the one loop 
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values at the above quoted Higgs masses, we will not discuss the values of 
these parameters in the MSSM. The reader should keep in mind that we 
do not know what the true form of the effective potential is near the phase 
transition, since we are interested in the case of a large Higgs mass for which 
the perturbative expansion breaks down. A more appropriate treatment may 
be the e-expansion as advocated in or some other technique based on 



renormalization group . However, it is not clear whether the e-expansion 
is accurate either in the case of the Standard Model or its supersymmetric 
version, so we do not know very well what the form of the effective potential 
is. The non-perturbative calculations [|12| , ||I5[ are still not at the level to be 
able to reconstruct the form of the potential for large Higgs masses. In order 
to illustrate our mechanism we made a simple assumption on the form of 
the effective potential, which is motivated by the one-loop finite temperature 
effective potential (with the cubic term, which gives rise to a first order 
phase transition, set to zero). This form of the effective potential should 
give a reasonable qualitative description the phase transition since it has the 
desirable generic features of a second order phase transition: it is maximally 
flat at the critical temperature where the quadratic term vanishes; the order 
parameter changes continuously as temperature drops. We believe that 
the use of the true form of the effective potential in the analysis that follows 
would not alter our main observation that Bq is suppressed by the expansion 
rate of the Universe at all times. 

The time dependence of C(t) and 9(t) in ( pO|) are both determined by 



4>{t), the value of the Higgs scalar expectation value. In Appendix B we 
show that to a very good approximation is given by the location of the 
minimum of the finite temperature effective potential (pi]): 



2D 

<PliT) = —[T^-T'] for T<To. (22) 

The relative departure of the true vev from this value, as is shown in Ap- 
pendix B, is suppressed by the expansion rate of the Universe. This remains 
true when the friction of the $ field is taken into account. 

In order to obtain the time dependence of 6, we make the simplified 
assumption 

^ - = e[l - TVToT^ (23) 
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where e is the net change in the CP violating phase (the relative phase 
between the Higgs doublets in the case of the two Higgs doublet model). 
Cline et al PU[ have done an extensive study of the dynamics of 6 and 



fields neglecting friction and find that for rather large range of parameters 
the above linear approximation is reasonable. As the analysis of Appendix 
B suggests, this conclusion should not be altered even when the friction of 
both fields is taken into account. 

In order to relate time derivatives to the derivatives with respect to tem- 
perature we will use the Friedmann equation for the radiation era: 

H=^ = h^,^, h^ = ^^, G = mpl = 106.25 (24) 

SO that tT^ = const and d/dt = —HTd/dT, or Hdt = —dT/T. 
^From Eqs. (^) and (p^ it then follows that 

-1/2 



1 T 

9 = He——— = eH— 

{To Ty - 1 To 



ToVT^-l , (25) 



and thus, maximizing C6, gives 



24 2D 
Bo = ^^^\gee\y^—HiTo)Tl (26) 



Using the equation 

27r2 

s = -^9*To' (27) 

for the entropy density, the baryon to entropy ratio is bounded from above 
by 

s - bV^TT^ 9* ' A To ^ ^ 

For To = llOGeV (M^ = 80GeV), the value of H/Tq is about 1.6 ■ 10"^^ 
and for To = 175GeV (M^ = MOGeV), H/Tq ^ 2.5 ■ 10-^^ ge ~ -1/4 
(we have assumed here that vi = V2), yt — l/"\/2- 2D/\ is bounded from 
above by about 7 when Mh = 80GeV. (For a larger value Mh = 140GeV, 
2D/X ~ 2.5.) Our final bound for the baryon to entropy ratio is then: 

— < 1.5 X lO^^^lel, (29) 
s 
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which completes the proof that without additional sources of out-of-thermal 
equilibrium such as topological defects, electroweak baryogenesis is not pos- 
sible if the electroweak phase transition is not first order. 

4 A lower bound on the baryon to entropy 
ratio 

We are now interested in obtaining an approximate lower bound on the 
baryon number today (which corresponds to -B = B{oo)). It proves con- 
venient to split the integral in Eq. (|T^ into two parts: the first for which the 
sphaleron rate F is large so that the time integral in the exponent is much 
larger than one, and the second part for which this integral is smaller than 
one. The time tn which separates the two regimes is defined by 

/ r{t")dt" = 1. (30) 

JtH 

On the first interval [tin, tn] "we can perform a partial integration to obtain 
B{oo) = ~Ce{tH)e^^ + Ce{t^n)e~ ^^^^'^'e"^ (31) 
+ r rft'^e-X'"^(*")^*"e-i - r dt'Cmme-I^'^^'"^''' (32) 

J tin at Jtjj 

We can immediately see that the second term vanishes if tm is chosen such 
that the Universe is initially in the symmetric phase. 

The idea behind the derivation of the lower bound is the following: the 
first and fourth terms have a negative sign, the third term is positive. We 
will show that the third term is smaller than the first in absolute value, and 
that hence the absolute value of the baryon number can be bounded from 
below by the absolute value of the final term. 

To find an upper bound on the third term, we replace the exponential 
factor by 1. The remaining integral exactly cancels the first term. Thus, the 
absolute value of B{oo) can be bounded from below by a lower bound on the 
absolute value of the last term. 

In order to bound the last term, it is necessary to know the exact tem- 
perature dependence of T (recall that for obtaining the upper bound on B{t) 
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this was not required). In the case of sphaleron-induced baryogenesis at the 
electroweak scale, we have: 



^(T) = 7T|/V^^ y 




4:71 



(33) 



where i?~ 1.5 — 2.7 is a slowly varying function of the coupling constant 
ratio, and 7 ~ 10~^fi;i, where ki is the one loop fluctuation determinant. 
Before we discuss the last term in (P2|), we will evaluate the integral in Eq. 



(|30D. We can solve this integral approximately by using the equation (|33|) for 
r and expressing the time differential in terms of the temperature differential 
with the help of the Friedmann equation (^41). The result is 



Xh 



r(xff) 

T{xh) 
T 



2D B Un H{To) 
A Xh V Oiw Tq 



(34) 



1 



a. 



B{gyX) V Att 



A /47r\- 



To 



2D \a^J B{gyX)H{To 



-X 



H 



which is solved for xh = 4>h/T ~ 1.2. This is nothing but the condition for 
the sphaleron erasure bound of Ref. |^ (see below). 

With the help of Eqs. (|13D, (||) and (gOj), the last term of Eq. can 
now be estimated as follows: 



tH 



dt'C{t')d{t')r{t')e 



■ IT nt")dt" 



> e~ 



tH 



dt'c{t')e{t')r{t')\{35) 



3 /2Dy/2 
2^ [~) 



\<^ge\yt 



Xh 



xu\/m + i] 



Bo 

e 



Comparing with Eq. it follows that the lower bound on the baryon to 
entropy ratio produced in the second order phase transition is only suppressed 
by a factor of e~^ compared to the upper bound. 



5 Conclusions 

In this paper we have shown that without additional sources which drive the 
system out of thermal equihbrium (such as a network of topological defects), 
baryogenesis during a second order electroweak phase transition is much too 
weak to be able to produce the observed baryon to entropy ratio. 
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We have shown that rather generically baryon number production in this 
case is suppressed by the expansion rate of the Universe. The case studied was 
that of a two Higgs doublet model with CP violation in the Higgs sector with 
the assumption that both Higgs fields couple the same way to the plasma and 
hence the phase transition occurs simultaneously for both Higgs fields. Even 
though in principle it would be interesting to study a more generic situation 
of a two stage phase transition, we beheve that the main conclusions of our 
letter still hold. 

In general, the source of CP violation may not be in the Yukawa terms, 
but in the neutralino and chargino mass matrices [^. In this case the 
effective shift in free energy should be proportional to a time derivative of 
the imaginary part of the masses, i. e. a gauge or Higgs coupling times a time 
derivative of the relative Higgs phase, so that the analysis is very similar to 
the one presented in this letter and again baryon production is suppressed 
by H. 

The second main result is that a non-vanishing baryon asymmetry is 
generated during such a phase transition. We have derived upper and lower 
bounds for the strength of this mechanism. They are of the same order of 
magnitude and suppressed by a factor H/Tq, where H is the Hubble expan- 
sion rate and Tq is the temperature of the phase transition. Note that the 
baryon production mechanism presented in this paper does not suffer from 
the sphaleron erasure problem since the bulk of baryons is produced when 
the sphaleron rate is switching off and becomes smaller than the expansion 
rate of the Universe. 

We have illustrated our mechanism for the effective potential which de- 
scribes a second order phase transition, such that the order parameter changes 
continuously at the phase transition but its (time) derivative is discontinu- 
ous, as opposed to crossover when the order parameter is a smooth function 
of time. However, from the point of view of baryon production in our mecha- 
nism, it is irrelevant whether the transition is second order or crossover since 
the bulk of baryon production occurs sufficiently below the phase transition 
when the difference between these two types of phase transition becomes 
immaterial. 

The formalism illustrated in this paper may be applicable to other baryo- 
genesis mechanisms taking place in a second order phase transition. Provided 
that a formula similar to Eq. ([I5|) holds, then an analysis similar to what 
was done in this paper would be applicable. The most optimistic value of 
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the resulting B/s is eH/To^/g^ ~ eTo/nipi, where as before e is a constant 
parameterizing the strength of CP violation in the particle physics sector de- 
termining the phase transition. Hence, provided that the scale of symmetry 
breaking is higher than about lO^GeV Q from the lower bound formula, it 
is possible to imagine a sufficiently powerful baryogenesis scenario involving 
only second order phase transitions. 



6 Appendix A: Derivation of the expression 
for the chemical potential 

In this Appendix we outline the derivation of Eq. ([T3| ) starting from Eq. (p!l|) 
in the limit ^ — > and rrii < T. 
Linearizing in 6, Eq. ([Tl|) gives 



where 



and 



^-PV^' + '"^ (37) 

1^1 + exp J x'^ + nij^Y 







x^dx expA/x^ + mf, 



x^ + (1 + exp Jx^ + 



^2 ^,2 



(38) 



The notation itlt = nii/T and x = p/T has been used. 

The integral can be performed explicitly in the limit 0, yielding 

C2 = 7r^/6. The second term is more tricky since the leading term (which is 
independent of uit) just represents a shift in /x. The method of evaluating 
I consists of breaking up the integration region into two intervals, the first 
from to where ^ y ^ 1, and the second from y to 00. In the first 
interval, we can set the exponential factors to 1 and evaluate the remaining 
integral to leading order in m-p, in the second interval we introduce a new 

''^A more realistic lower bound on the scale of symmetry breaking is about 5 x 10^°GeV; 
this is easily obtained from Eq. (p9|) and the requirement e < 1. 
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integration variable z = ya;^ + mj,, obtaining 

^JW^T ^^(exp;z+ 1)2' ^ ' 

The second integral (coefficient m|.) can be performed explicitly and is of 
the order mf^, the ffist can be performed by ffist integrating from to cxo 



and then subtracting the result obtained by integrating from to y?/^ _|_ 772^,. 
The ffist can be done explicitly, and gives a constant which can be viewed 
as a shift in /x, the second integral can be approximated to leading order 



in \Jy'^ + rrij^, and to this order exactly cancels the leading y dependence of 
integration of I from to y. A good computational check is that the final 
result must be independent on the cut-off y, since the cut-off is introduced 
solely for computational convenience. Thus, to leading order in mx, the 
result for / is 

/ = C2-^(1-^). (40) 

The first term in J should have no physical effect since in the zero mass 
limit Q is pure gauge and must have no physical effect whatsoever. It can be 
viewed as a shift in the chemical potential [i: 

/i = /i =F geQ, (41) 

where Jl denotes the physical chemical potential with a correct physical in- 
terpretation: /2 is proportional to the particle density perturbation. With 
the substitution of Eq. (|1T|), and inserting the above results for J and J^, 



Eq. (p6| ) yields the result for [i given in Eq. (pTsp. (Note that in the main 
text we have denoted the physical (shifted) chemical potential by /i.) 



7 Appendix B: The irrelevance of dissipative 
effects 

In this Appendix we will show that in the case of a second order phase 
transition dissipation does not effect the evolution of the field and it is 
irrelevant for baryogenesis considerations. This is in contrast to the case of a 
first order phase transition where it is known that dissipation is crucial 
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for the dynamics of the transition and baryogenesis. We will study the case 
of the (p field, but an analogous analysis applies to the relative Higgs phase 

e. 

The equation of motion for the Higgs condensate is 

4> + 3H^ + ^^^+^^-^^ = 0. (42) 

The form of the dissipative term 7,^0 can be obtained in the WKB approxi- 
mation using the method developed in [^. Assuming that the main source 



of dissipation comes from the tree-level process in which the Higgs particle 
scatters off the top and decays, we obtain for the friction 

where F,^ ~ 0.3a^|/^T is the Higgs decay rate. The lessons to learn from 
this simple calculation are the following. Firstly, the friction term is local 
in time. This is to be expected since the time scale on which the Universe 
evolves ~ 1/H is huge in comparison to the time scale on which perturbations 
are damped r^p ~ l/F^. This is in contrast to what happens at a first order 
phase transition where the system is perturbed out of equilibrium in a non- 
local manner in the sense that perturbations are sourced and then propagate 
across the phase boundary. Secondly, in the case of one fluid, friction is 
related to viscosity, while in our case of many fluids interacting, the main 
source of friction are couplings between different species. Thirdly, the friction 
for the field is immense in comparison to the expansion rate of the Universe. 

Since the time scale on which the Universe evolves is given by ~ 1/H, 
we can neglect the first two terms in Eq. (|42|) . We now write Eq. (|42|) in the 
following form: 

= <j)o + ^, 7^0 + 8D(r2-T2)y, + 3[2DA(r2-r2)]iV + Ay.=' = (44) 
The solution is 

_ -f^H 11 

f - -'^0^^(T2/T2-1)2 ^^^^ 

which means that the instantaneous value of is very well approximated 
by the equilibrium value 0o in the sense that the relative correction ip/(f) is 
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suppressed by the expansion rate of the Universe. One can check that indeed 
all of the following 0, ip^ and ip^ are suppressed by at least H in comparison 
to which justifies the above approximate solution. 

We have thus shown that friction does not alter in any significant manner 
the evolution of the Higgs condensate and the use of 0o in the main text was 
justified. An analogous consideration applies to the 6 field. 
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